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ABSTRACT 

We  investigate  the  stability  of  two  phase  Couette  flow  of  different  liquids  bounded  between 
plane  parallel  plates.  One  of  the  plates  has  a  time  dependent  velocity  in  its  own  plane,  which 
is  composed  of  a  constant  steady  part  and  a  time  harmonic  component.  In  the  absence  of  time 
harmonic  modulations  the  flow  can  be  unstable  to  an  interfacial  instability  if  the  viscosities 
are  different  and  the  more  viscous  fluid  occupies  the  thinner  of  the  two  layers.  Using  Floquet 
theory,  we  shov;  analytically  in  the  limit  of  long  waves,  that  time  periodic  modulations  in 
the  basic  flow  can  have  a  significant  influence  on  flow  stability.  In  particular,  flows  which 
are  otherwise  unstable  for  extensive  ranges  of  viscosity  ratios,  can  be  stabilized  completely 
by  the  inclusion  of  background  modulations,  a  finding  that  can  have  useful  consequences  in 
many  practical  applications. 
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1  Introduction 


It  is  well  known  that  plane  Couette  flow  of  two  superposed  fluids  of  diflerent  viscosities 
can  be  unstable.  This  instability  is  absent  when  the  viscosities  are  equal  and  is  a  result  of 
an  interfacial  deflection  growing  due  to  viscosity  stratification.  This  niechanisin  was  flrst 
described  by  Yih  (1967),  who  used  perturbation  theory  to  obtain  an  analytic  expression  for 
the  grow’th  rate  in  the  long  wavelength  limit.  The  results  can  be  suniniarized  as  follows  (for 
clarity,  in  the  absence  of  gravity):  Instaoility  is  possible  only  if  the  viscosities  are  different. 
Arrangements  where  the  more  viscous  fluid  occupies  a  thinner  layer  than  the  less  viscous  one 
are  unstable,  whereas  converse  arrangements  are  linearly  stable  to  long  waves.  Tlie  insi  ability 
is  present  at  all  Reynolds  numbers,  although  in  the  limit  of  zero  Reynolds  number  the  growth 
rates  become  asymptotically  small.  Flow,  therefore,  is  necessary  to  excite  growing  waves. 
Moderate  surface  tension  is  negligible  to  leading  order  for  long  waves,  but  arts  to  stabilize 
very  short  waves.  Yih’s  work  was  extended  more  recently  by  other  investigators.  Hooper 
and  Boyd  (1983)  considered  the  linear  stability  of  two  co-flowing  viscous  liquids  of  different 
viscosities  separated  by  an  interface  and  extending  to  infinity.  This  is  Yih's  problem  in  the 
absence  of  walls  and  can  be  used  as  a  useful  model  in  the  study  of  short  waves  (short  waves 
at  the  interface  would  not  feel  the  effect  of  the  walls,  to  leading  order  at  least).  Instability 
is  found  in  ail  cases.  Long  waves  are  stable  for  ail  viscosity  ratios,  while  .short  waves,  in 
the  absence  of  surface  tension  are  unstable  with  asymptotically  small  growth  rates  -  when 
included,  surface  tension  damps  out  all  waves  which  are  short  enough.  Disturbances  with 
general  wavelengths  are  unstable  and  a  maximum  growth  rate  is  attained  at  some  viscosity 
ratio.  This  work  was  extended  by  Hooper  (1985)  who  included  a  bounding  moving  wall  for 
the  lower  fluid.  The  flow  is  shown  to  be  unstable  to  long  waves  as  long  as  the  bounded  layer 
is  also  the  more  viscous  one,  in  agreement  with  Yih’s  findings.  An  explicit  expression  for 
the  eigenvalues  is  also  given  in  terms  of  Airy  functions. 

The  stability  of  Yih’s  problem  for  general  disturbances  and  including  density  differences 
and  surface  tension,  was  first  solved  numerically  by  Renardy  (1985).  She  presents  specific 
examples  that  show  that  the  general  long  wavelength  result  described  by  Yih  and  Hooper 
and  Boyd,  can  be  extended  to  interfacial  waves  of  arbitrary  wavelengths.  Another  interesting 
feature  is  the  finding  of  a  second  mode  of  instability  at  relatively  high  Reynolds  numbers; 
this  second  mode  has  long  wavelength  and  growth  rates  comparable  to  the  first  interfaciaJ 
mode.  Renardy  (1987)  carried  out  an  analytical  study  of  the  s^me  problem  but  for  fluids 
with  slightly  different  mechanical  properties.  Analytical  expressions  for  the  eigenvalues  are 
constructed  and  instability  is  established  in  the  thin  layer  limit  if  the  thin  layer  is  occupied 
by  the  more  viscous  fluid. 
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The  stability  ohararteristirs  of  oscillatory  Hows  are  iiiuch  more  r<;jiipli(  ated  than  those 
of  their  steady  counterparts.  There  are  many  examples  wliere  the  inclusion  tT  an  os<  i!lalorv 
component  to  the  steady  flow  can  enhance  or  reduce  stability  (for  plane  I’oiKeuille  flow,  fur 
example,  see  Crosch  and  Salwen  (1968),  Hall  (1975),  von  Kerrzek  (1982);  fur  a  review  on  the 
stability  of  oscillatory  flows  s^e  Davis  (1975)).  As  far  as  we  know  the  stability  of  oscillatory 
two-phase  viscous  flow  has  not  been  studied.  Yih  (1968)  considered  the  stability  of  a  viscous 
fluid  layer  on  a  flat  plate  performing  a  simple  harmonic  motion  in  its  own  plane,  Tlnue  is  no 
viscosity  stratification  (the  upper  fluid  is  air)  and  the  flow  in  the  absenct*  of  the  oscillation 
contains  no  steady  velocity  component  and  is  linearly  stable.  Yih  showed,  using  a  long  wave 
Fioquet  theory,  that  the  oscillatory  How  can  become  unstable  even  though  in  the  ab.sence 
of  oscillations  the  flow  is  stable.  This  was  extended  by  von  Kerczek  (1987)  to  flow  down 
a  vertical  plate  which  is  performing  a  simple  harmonic  motion  in  its  own  plane.  This  flow 
is  unstable  even  in  the  absence  of  oscillations  and  von  Kerczek  uses  Yih’s  long  wavelengtii 
expansion  to  establish  windows  of  instability.  Our  interest  is  to  apply  si]ch  a  study  to  an 
oscillatory  two-phase  flow.  In  the  problem  we  study  here,  depending  on  tin*  fluid  occupation 
areas  in  the  unperturbed  state  (see  earlier  comments),  the  flow  can  be  stable  or  unstable.  In 
the  former  case,  therefore,  a  quasi-static  approach  would  yield  background  velocity  profiles 
for  each  instant  in  time  (time  is  treated  as  a  parameter  in  qua^i-steady  theories)  which  are 
linearly  stable  while  in  the  latter  instance  instability  ensues  for  all  parametrized  profiles.  A 
more  yielding  approach  is  that  of  Fioquet  theory  where  stability  or  instability  is  judged  on 
the  overall  growth  or  decay  with  time  of  a  perturbation  over  a  complete  period  of  the  forced 
oscillation.  Analytically  this  means  that  time-periodic  solutions  are  constructed  which  can 
amplify,  remain  neutral  or  decay  depending  on  whether  the  Fioquet  ex])onent  is  i)ositive. 
zero  or  negative.  We  present  representative  results  for  several  cases  which  in  the  absence 
of  background  oscillations  the  flow  is  unstable  while  imposed  oscillations  can  stal)i!ize  the 
flow  completely.  On  the  other  hand,  flows  which  are  stable  can  become  unstable.  We  note 
that  all  the  results  given  here  are  in  terms  of  explicit  but  long  fornnilae  wliich  yield  stability 
results  with  little  computational  effort.  Such  explicit  solutions  are  of  considerable  value  in 
the  generalization  of  the  stability  problem  to  arbitrary  wavenutnbers  by  use  of  continuation 
methods  for  instance. 

The  article  is  organized  as  follows.  Section  2  derives  the  undisturbed  flow  as  ar  I'xact 
shear  flow  .solution  of  the  Navier-Stokes  equations.  In  Section  d  the  liin-ar  stability  problem  is 
formulated  and  the  interfacia!  conditions  are  written  down.  Explii  it  solutions  ar('  dcweloped 
by  an  expansion  procf 'hire  carried  out  to  three  orders,  so  that  we  may  d('t('rmine  the  first 
stage  where  the  expansioii  of  the  Flocpiet  exponent  is  a  o'al  tenn  In  S--  tion  1  we  prese’,, 
the  results  of  the  st...bility  analysis  and  finally  in  .St'ction  5  we  flraw  some  conclusions. 


2  The  basic  state 


Two  incompressible  fluids  of  equal  densities,  p,  and  different  viscosities,  occupy  a  region 
of  depth  L  between  parallel  infinite  flat  plates.  The  lower  plate  is  fixed  aijd  the  upper 
plate  moves  in  its  plane,  with  a  steady  velocity  I'o,  together  with  superposed  sinusoidal 
oscillations,  so  that  its  velocity  at  time  t”  is  given  by  f  ^o  +  A  cos  (u.*/’).  The  fluids  ar<‘  assiinieil 
to  be  immiscible  and  form  separate  layers,  the  upper  fluid  has  viscosity  pi,  tin*  lower  fluid 
has  viscosity  p.^  and  a  depth  D.  We  denote  the.se  as  regions  (/)  and  {II)  respectively 
We  look  for  exact  solutions  of  the  Naviei-Stokes  equations  for  a  constant  interface  position 
and  velocity  component  in  the  streamwise  direction  alone  which  depends  on  time  and  the 
vertical  coordinate  only.  Using  cartesian  coordinates  the  exact  flow  is  described  l>y 

the  following  equations 


dWl  p,  dWi 

dt*  p  dy*^  ' 

dul  _  P2 

dt*  p  dy*^  ’ 


L>y->D  (/), 

D>y‘>0  {II). 


The  boundary  conditions  are,  no  slip  at  the  walls, 


U]{y*  =  L)  =  Uo  + A  cos  (ujt*), 

TTAy-^O)  =  0. 


and  continuity  of  velocities  and  stresses  at  the  interface, 


IT.iy-^D)  =  TT^iy^D), 

=0)  =  =0)- 

The  solution  is  easily  found  by  separation  into  a  steady  part  and  an  oscillatory  part  of  ihe 
same  frequency  as  the  plate  motion.  The  basic  flow  is. 


IT, 


f^pUoy*  - UoD{p^-^)  ^  ^ 

Lp2-  D{P2- li\)  ^  ^ 


Lpi  -  D{p2-  fii) 


+  A7Z  <  2K  sinh 


0‘y' 


T713 


where  71  denotes  the  real  part,  (i*  =  {ipijjfp\)^  and  we  define  the  ratio  of  the  lower  fluid 
viscosity  to  that  of  the  upper  layer  by  the  parameter  m  =  P2/pi-  Tfi‘'  constants  A’,  bj  and 
Li  are  given  by 


K  - 


1 


2  sinh  {  co.sh(/f* 


{D- 


L))  +  2/n2  cosh 


sinh  (ii*  (fl—  L)) 


i) 


Lr 


A'exp(-/i*D) 


rn^  cosh 


+  sinh 


m2  / 


(2) 


Li  = 


A^exp(/rD) 


m2  cosh 


(3) 


The  flow  is  non-diniensionalized  next  with  D  as  lengthscale  and  Uq  as  velocity  scale,  while 
time  r  and  pressure  P'  are  non-dimensionalized  by  DjUo  and  plJ^  respectively.  This  intro¬ 
duces  the  following  non-dimensional  groups  (used  later  also); 


A  =  A,  ^  =  ,3-D, 

pDUi'  Co' 

The  parameter  a  measures  the  relative  depths  of  the  two  fluids,  A  measures  the  magnitude 
of  the  externally  imposed  oscillations,  /3  is  a  measure  of  the  Stokes  layer  thickness  due  to 
the  oscillatory  flow,  Re  is  the  Reynolds  number,  and  Q  is  the  non-dimensional  frequency  of 
the  plate  oscillations.  Under  these  scalings  the  non-dimensional  no-slip  condition  becomes 


a  = 


Re  = 


L 

D’ 

UpDp 

/^i 


(/i  (y  =  a)  =  1 -f  Acos(flt) , 

F2(y  =  0)  =  0, 


and  the  basic  flow  is: 


Ui 


my  —  m  -h  1 
am  —  m  -f 

y 

am  —  m  -f  I 


-  -h  A7^{(Tle(^^)  -b  , 


-f  A7?.  |2/('  sinh 


0y 


,  m2 


o(>nt) 


The  constants  K ,  L\  and  Li  appearing  above  are  in  their  non-dimensional  foims  which  are 
readily  available  from  (l)-(3)  by  replacing  (PD  by  (L  When  A  =  0  we  recover  the  basic 
(louette  flow  of  different  fluids  in  a  channel  (see  Yih  (1967)). 


3  Stability  equations. 

The  equation  governing  the  linear  stability  of  parallel  shear  flows  is  the  Orr-Sommerfeld 
equation  and  arises  from  linearization  of  the  Navier-Stokes  equations  and  a  norma!  tiiode 
expansion.  In  the  present  problem  the  background  flow  is  time-dept'udent  and  so  tln^  eigmi- 
functions  depeinl  both  on  t  and  y.  In  particular,  the  pertnrbaiion  sireanifnnct ion  is  tak(>n  to 
have  the  form  <?h,^  {•P-y.  i)  —  <l>]a  where  o  is  the  wavennmlx'r  of  the  di.sturl)ance  and 

sul«scripts  1,2  denote  quantities  In  regions  (I)  and  (II)  respectively.  Th('  stability  eciiialions 
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in  regions  (I)  and  {II)  are,  then, 


'd  ,  —  \( ■2\  ‘  ■  1  ^ 


^Vl  „  ,4 

+  a  '?! 

O  ,  4 

■  2a^ -^T-r  +  «  <P2 
dy^ 


dy* 

d%2 

dy* 


(4a) 

(4b) 


The  linearized  boundary  conditions  are 


<l>i  (a)  =  0=  ^{c). 

^2(0)  =  0=^(0). 


(5) 

(6) 


We  now  let  the  interface  position  be  y  =  1  -t-  Sh{t)  6*°^  where  S  «  1  is  the  infinitesimally 
small  amplitude  of  the  perturbations.  The  linearized  interfacial  conditions  are  found  by 
expansion  of  the  exact  conditions  about  y  =  1  in  powers  of  6  and  retention  of  the  leading 
order  contributions: 


<f>\  =  d*2» 
dy 


d(i>2  .  (I  -m)dUi 

dy  m  dy  ' 


p  +  =  m 


dy 

-  R, 


^^2  ,  2 


+  Q!>2 


dy^ 

'd^2 


dy 

^d(i>2 


dy 


~  ^  \  TT  ~ 

\  dy^  dy  dy  ^ 


-  Re 


[dy^ 

d^r  ,  .  rjdh  .  ,  dUr 

dydt  dy  dy 

0^2  ,■  77  d<^2  .  .  dUr 


+ 


ia^/17 

"X"' 


The  linearized  kinematic  condition  at  y  =  1  reads 

— — f-  iahU\  +  iot<f>\  =  0. 

at 


(7) 

(8) 

(9) 


(10) 


(ii; 


Equations  (4-11)  constitute  the  partial  differential  system  together  with  boundary  condi¬ 
tions  that  governs  the  stability  of  the  flow.  The  problem  contains  time  and  one  space  variable 
and  in  general  requires  integration  in  time  as  an  initial  boundary  value  problem.  Stability  or 
instability  is  determined  by  the  large  time  evolution  of  initial  perturbations  corresponding 
to  different  wavenumbers. 

In  this  work  we  proceed  analytically  by  studying  the  stability  characteristics  of  long 
waves  (a  <<  1).  Such  analysis  provides  a  significant  amount  of  information  since  it  enables 


development  of  closed  form  expressions  which  can  be  used  to  check  whether  unsteadiness  in 
the  background  flow  can  affect  the  evolution  of  linear  perturbations.  In  the  ('ouette  j^roblem 
studied  by  Yih  (1967),  instability  was  established  for  all  Reynolds  numlx'rs  by  an  analogous 
long  wave  expansion.  The  analysis  of  long  waves  is  also  useful  in  guiding  numerical  solutions 
of  the  initial  value  problem  at  general  wavenumbers. 


3.1  Solution  for  long  waves 

We  now  proceed  analytically  by  considering  the  asymptotic  limit  a  0  which  corresponds 
to  disturbances  with  long  streamwise  wavelengths.  We  use  a  method  simila.r  to  the  Floquet 
theory  applied  by  Yih  (1968)  to  obtain  analytically  the  growth-rate  of  the  disturbance  at 
the  interface.  We  look  for  solutions  for  the  streamfunction  and  interfacial  amplitude  as  a 
power  series  in  a  1  of  the  form,  for  j  =  1,2, 

<f>]  (2/)0  =  (j/iO  +  +  •  •  •)exp((^o  +  a'6'i  +  . .  .)t)exp{ikat.) , 

h{t)  =  {ho{t)  +  ah^  (t)  +  . . .)  exp[{6Q  +  a6]  .  .)t)exp  (ikat) , 

k  ~  ^0  T  Oik]  -f- . .  • 


where  ko,  k^,. . .,  9q,0],.  . .  are  real  constant®,  and  /lo  (<) ,  (f ) , . . .  are  periodic  in  time  t. 

By  writing  the  streamfunction  and  interfacial  deflection  in  this  form,  the  linear  long  wave 
stability  of  the  unsteady  flow  can  then  be  determined  by  calculating  the  first  non-zero  9. 
which  corresponds  to  exponential  growth  or  decay  of  the  disturbance  to  the  basic  flow. 

We  now  consider  the  system  of  governing  equations  (4a-b)  and  the  boundary  and  interface 
conditions  (5-11)  to  leading  order,  0(1)-  Firstly  the  kinematic  condition,  (11)  yields 

^  =  (12) 

Now  since  9o  is  real  and  ho  is  periodic,  we  have  two  choices;  we  may  take  fio  =  0  so  that  9o 
is  as  yet  unprescribed,  or  alternatively  0o  =  0  and  then,  without  loss  of  generality,  ho  —  1. 
We  shall  in  fact  assume  the  latter,  since  ho  —  0  leads  to  damped  modes  only  as  shown  later 
in  Section  3.2. 

After  making  the  transformation  i^jo  —  <l>jo  +  for  each  region,  j  —  1,2.  the  leading 
order  system  becomes 


1 

dy'^dt 

'  R~e 

dy^  ’ 

m 

dy'^dt 

w. 

dy^ 

-  aj)  = 

--  ‘  +  7F‘‘' 

(13a) 


(13b) 


(bla) 
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^^20  {y  =  0,  It)  =  0, 

dy  "am-m+1 


+  [^  (Z„e^“-L,e-^»)  e*«‘  +  ..cl 


oy  am  —  m -j- 1 

Here,  c.c.  denotes  the  complex  conjugate, 

At  the  interface  t/  =  1  we  have 


(14b) 

(14c) 

(14d) 


010 

=  020, 

(1.5a) 

C?01O 

dli>20 

(15b) 

dy 

dy  ’ 

5^010 

5^020 

(15c) 

dy^ 

=  “  V  ' 

1 

<9^010 

dy^ 

p  ^Vio 
'  dydt 

II 

3 

0 

D  <9^020 
'  dydt  ’ 

(1.5d) 

These  equations  admit  solutions  of  the  form: 

^10  =  Aio  (y  —  a)^  4- J?io  (y  ~  a)^  + - - - V  ^ 

am  —  m  +  1 

+  {[Cio  sinh  (/3y)  +  D^o  cosh  (/?y)  +  Eioy  +  Fio]  e’”‘  +  c.c.}  , 

020  =  A2oy^  +  B^oy^  + - - - - 

om  —  m  +  1 

+  { [C2osinh  (/Jj/rn'i}  +  D2ocosh  (^j3ym~i^  E-zoy  +  Fzo  e’^‘  4-  c.c.}  . 

The  twelve  constants  are  found  analytically  by  substitution  into  the  above  boundary  and 
interface  conditions,  the  resulting  equations  which  determine  Ajo,  A20,  Fjo,  F-zq  are  given 
in  the  appendix. 

At  the  next  order,  0(a),  the  kinematic  condition  (11)  becomes 

dh\ 

-^  +  01  +  1010  (l,  0  +  =  0  (16) 

Our  solution  strategy  in  obtaining  the  eigenvalues  analytically,  is  based  in  part  on  the  per¬ 
turbation  eigenfunctions  and  interfacial  amplitude  being  time  periodic.  In  equation  (16) 
above,  0io(l,t)  is  the  sum  of  a  real  constant  and  a  time  periodic  function,  while  A-j ,  and 
01  are  real.  In  order  to  obtain  periodic  solutions  hi{t)  the  following  constraints  need  to  be 
satisfied, 


■010  (U  ~  ^10  (c-  —  1 )'  ~  B\o  (1  —  n)^  -f 


m  («  —  1 ) 
am  —  m  4-  1 


{ 


^  [Ciosinh  (0)  +  Dio  cosh  {0)  +  Dio  +  Dio] 


C.C., 


(17) 


where  the  superscripts  (s)  and  (us)  denote  the  steady  and  unsteady  parts  respectively. 

Before  solving  the  0(0;)  system,  let  us  first  consider  the  O  {a^)  kinematic  condition, 
which  can  be  written  as 


dh,2 

dt 


9-2  ikih\  +  ik'i  +  thxU \  +  2<^i]  ( 1 ,  f )  —  0. 


Using  X  to  denote  the  imaginary  part,  we  see  that  for  h2{t)  to  be  periodic  in  t  we  require, 

0,  =  j|[/iiFi+<Aii(1)]^'^}. 

Ulearly  then  it  is  sufficient  to  solve  for  the  steady  part  of  the  eigenfunction  (p^x  to  0{o:) 
in  order  to  determine  62  and  hence  the  linear  stability  of  the  interface.  We  note  that  at 
this  order,  products  of  functions  which  have  a  time  dependence  given  by  exp(±iDt),  yield 
additional  steady  terms  in  the  governing  equations  and  boundary /interface  conditions.  The 
eigenvalue  82  is  therefore  determined  by  both  the  steady  and  time  oscillatory  motion  of  the 
interface.  The  time  independent  perturbation  equations  are  given  by 

i(d 


1 

De  dy^ 


{^’1  + 


K  dy^ 

They  admit  a  particular  solution 


d^xo  dWx 
dy'^  dy^ 

0^2 


<Pic 


(1.  ,  T7\^  ^  ^ 

i*'  +  -w " 


I  I  ^ - *»1 


(») 


{18a) 

(18b) 


dy 


X21  = 


dy 


dy 

- 


(s) 


TF  ^  ^V^20  ,  dU2 


dy 


dydydy, 

(d 

dydydy. 


(19a) 

{19b) 


so  that  the  general  solution  of  equations  (ISa-b)  is 


<Px\  ~  'dll  (j/ —  1)^  +  Du  (j/ —  1)'^  +  Dn  (?/ —  1)  +  Dll  +  \ii, 

4>2\  ~  'd-ii  (y  —  1  )^  +  D21  iy  —  1  )'^  +  (-21  iy  —  1 )  +  On  +  \2!  • 

The  eight  constants  are  found  by  applying  the  following  boundary  and  interface'  conditions: 

4'’i’  (?y  =  «)  -  0, 
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(y  =  0) 


<l>n 

941 

dy 

dy^ 

dy^ 


(»)  _ 


0, 

0, 

0, 


(rn-  1) 


=  m 


—  m 


dy 

ay2  ’ 

9'Ul1 

dy^ 


771 


h, 


dUx 


dy 


i>) 


-R. 


ikid4>\Q  .jj  d(f>\Q  .  j  dU\ 


dy 


+  Rt 


dy 


dy  J 


(») 


ikid(i>2Q  ,  -TT  ^^20  .  , 


The  above  solutions  and  boundary  conditions  provide  a  system  of  eight  inhomogeneous 
algebraic  equations  for  the  eight  constants  Axi,A2\,---,Dxi,D2\  and  are  solved  explicitly 
here,  (see  appendix). 


3.2  Damped  disturbances 


Before  presenting  the  results  of  the  calculation  of  the  eigenvalue  0^,  we  show  that  the  possi¬ 
bility  ho  =  0  yields  damped  waves.  To  achieve  this  we  write  the  leading  order  streamfunction 
as 


<f>jo  =  Qj  iy)  exp  (^of ) ,  where  j  =  1 , 2. 


This  satisfies  the  equations 


RgOo 

Re^O 


dy"^ 

<PQ2 

dy^ 


tRl 

dy* 

d*Q2 

m-r-r 


which  are  solved  subject  to 

13.(0)  =  ^(0)  =  0, 
0j(»)  =  ^(a)  =  0, 

(3  (1)  =  0!(1), 


(20) 

(21) 
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<ly 


iy 


m— r-r-  ( 1 ) . 


dy^  dy^ 

dy 


^(l)-ft,0o^(l)  =  in-^{\]-HA^ 

•J  n'^  nt!^ 


dy^ 


dj^ 

dy 


(!)■ 


Multiplying  (20)  and  (21)  by  ^  and  ^  respt^ctively,  integrating  between  0  and  a,  and 
imposing  the  above  boundary  and  interface  conditions  yields 


fif&o  ~  ~^n 


fo  (^)  +  /r  (^)  dy 


[/o(^^)  ‘^y  +  /r(^)  dy 

Hence  all  disturbances  are  damped  for  this  case. 


<  0. 


4  Results 

The  non-dimensional  parameter  A  is  the  magnitude  of  the  sinusoidal  oscillations  of  the 
upper  plate  relative  to  the  magnitude  of  the  steady  velocity.  When  A  =  0,  this  problem 
reduces  to  that  of  steady  plane  Couette  flow  of  two  superposed  fluids  of  differing  viscosity, 
(as  characterized  by  the  ratio  m). 

As  an  initial  check  on  the  validity  of  our  analysis  we  firstly  consider  this  steady  case,  which 
corresponds  to  the  problem  solved  by  Yih  (1967).  In  order  to  make  a  direct  comparison  of  the 
results  we  must  first  note  that  the  length  scale  used  by  Yih  (1967)  is  equivalent  to  (L  —  D) 
whereas  in  this  work  we  use  the  lower  fluid  depth  D.  Yih  also  defines  a  parameter  n  =  d^/di 
the  depth  ratio  of  the  lower  to  upper  regions,  a  brief  calculation  shows  that  u  =  (o  —  1 )"’ .  In 
view  of  these  notational  differences  the  Reynolds  number  and  streairiwise  wavenumber  are  not 
the  same  as  their  counterparts  defined  in  this  work,  we  must  first  make  the  transformations 
ctyih  — >  na  and  Ryiu  — »  n/?^,  we  hence  write  J  (m,a,  A  =  0)  =  {a  -  1  )~^  02- 

(Calculation  of  and  hence  J,  is  in  theory  an  analytical  task,  in  practice  however,  the 
twenty  constants  obtained  in  the  integration  of  the  0(1)  and  0(rt)  systems,  are  lengthy 
expressions  which  are  found  most  efficiently  by  u.se  of  a  symbolic  aigel>ra  package'.  PYr 
completeness,  the  simultaneous  equations  which  determine  these  constants  are  given  in  the 
appendix. 

For  a  given  visrrxsily  ratio  r/t,  depth  ratio  a.  frerpiencv  fl  and  magnitude  A.  the  value  of  O2 
is  calculated.  When  the  real  Flotpiet  exponent  Oi  >  0  the  disturbance  to  tlu'  lia.sic  flow  wifi 
grow  exponentially  in  time  ami  th<’  flovv  will  be  unstalile,  .similarly  negati\’e  02  corre.spomls 
to  a  damped  interfacial  disturbance,  and  the  flow  is  said  to  be  liiu'arly  stable. 
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Figures  1(a)  and  1(b)  plot  the  values  of  J  against  m  for  a  steady  flow,  (A  —  0),  when  tlie 
lower  fluid  is  more  viscous  than  the  upper  layer.  We  have  chosen  equivalent  depth  ratios  to 
those  illustrated  by  Yih  (1967),  the  results  are  identical. 

Figure  i(a)  corresponds  to  a  flow  in  whicii  the  uppei  layer  i.s  deeper,  for  a  mure  viscous 
lower  layer  the  growth  rate  is  positive  and  the  flow  is  always  unstable.  With  a  shallower 
upper  layer,  Figure  1(b)  indicates  that  there  is  a  region  of  stability,  depending  on  the  size  of 
the  viscosity  ratio,  although  for  sufficiently  large  in  the  flow  will  become  unstable. 

For  unequal  depths  of  fluids  the  growth  rate  tends  linearly  to  zero  as  m  approaclies  unity, 
since  this  is  a  hidden  mode  for  flows  of  equal  viscosities. 

To  obtain  results  for  a  less  viscous  lower  region,  we  have  verified  the  following  transfor¬ 
mation  given  by  Yih 

mj  (m,  n,  A  =  0)  =  n^J  f  ,  A  =  o')  ,  where  n  =  (a  —  1  )“V  (22) 

\m  n  J 

The  results  for  the  steady  flow  when  the  more  viscous  fluid  occupies  the  upper  region  can 
then  be  inferred  from  Figures  1(a)  and  1(b). 

Before  discussing  the  results  for  the  full  time-oscillatory  problem,  let  us  first  investigate 
this  property  in  more  detail.  Yih  observes  that  the  phase  speed  of  the  disturbed  flow, 
(cq  using  liis  notation),  is  equal  in  magnitude  and  opposite  in  sign  vvlien  the  viscosities 
and  depths  of  the  fluid  layers  are  interchanged.  This  is  because,  the  original  flow  may 
be  recovered  identically  by  making  a  Galilean  transformation  from  the  coordinates  {x,y,t) 
to  where  ^  =  x  -  t,  hence  (|^,  (^,  J  -  The  upper  boundary 

now  becomes  stationary,  (relative  to  the  moving  frame  of  reference),  and  the  lower  plate 
moves  with  constant  velocity  —1.  Since  gravity  is  neglected  for  this  problem,  inverting  the 
geometry  leaves  the  eigenvalue,  and  hence  the  stability  characteristics,  unchanged.  In  order 
to  compare  the  growih-rates  quantitatively  however  we  must  note  tliat  now  the  Reynolds 
number,  Re  =  UoDp/fii,  and  lengthscale  D  must  be  rescaled  on  the  neir  upper  fluid,  hence 
the  multiplying  factor  n^lm,  as  given  by  equation  (22)  above. 

We  now  consider  the  problem  when  A  7^  0,  the  flow  is  dependent  upon  the  time  t  and 
analogous  conclusions  from  a  change  of  reference  ficone  do  not  hold  and  so  for  completeness 
we  also  need  to  consid.T  arrangements  when  when  the  lower  fluid  is  less  viscous  than  the 
upper  layer.  Since  the  unsteady  contribution  is  periodic,  we  need  only  consider  A  >  0.  In 
fact  it  is  clearly  seen  that  A  is  a  multiplicative  factor  in  the  basic  flow,  O10.O20  and  also  h^. 
and  these  terms  contribute  to  O2  as  products,  thus  —A  gives  identical  results. 

When  the  lower  fluid  is  more  viscous,  rn  >  1,  the  interface  is  destabilized  by  the  intro¬ 
duction  of  an  unsteady  basic  flow.  Figure  2(a)  shows  that  for  a  =  .'ho  the  already  unstable 
steady  mode,  (broken  line),  is  made  more  unstable  upon  increasing  the  magnitudt'  of  the 
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oscillations.  Similarly,  when  the  lower  fluid  is  deeper,  for  example  a  ~  1.4,  s(j  that  in  the 
absence  of  background  time  periodic  modulations  the  flow  is  stable  for  moderate  in  at  least, 
the  interface  is  again  destabilized  as  indicated  by  Figure  ‘2(b).  The  value  of  m  below  which 
the  flow  is  stable,  decreases  as  A  increases,  for  example  when  A  =  0.2.  the  neutral  distur¬ 
bance  is  obtained  for  a  viscosity  ratio  m  m  86.2319,  wdiereas  for  A  =  0.4,  m  2:  66.2506.  This 
is  shown  more  clearly  by  the  neutral  curve  in  Figure  3:  if  the  magnitude  of  the  oscillations 
increases  beyond  the  critical  value  A  2r  0.5695,  the  interface  is  unstable  for  all  m  >  1. 

When  the  lower  fluid  is  less  viscous,  in  <  1,  the  results  are  more  signiflraut.  For  all 
depth  ratios  we  find  that  the  time-dependent  oscillatiojis  dramatically  stabilize  the  interfacial 
disturbance.  For  a  deeper  upper  layer,  for  example  a  —  2.25  as  in  Figure  4(a),  the  range  of 
viscosity  ratios  m,  for  which  the  real  growth  rate  is  negative  is  increased  for  the  unsteady 
background  flow  and  as  A  is  increased  the  flow  is  stabilized  further.  Figure  4(b).  illustrates 
how’  the  unstable  steady  mode  corresponding  to  a  shallow  upper  layer,  can  be  completely 
stabilized  provided  A  is  made  large  enough.  Such  behavior  is  shown  collectively  in  Figure  5 
which  can  be  used  to  predict  oscillation  amplitudes  wdiich  completely  stabilize  the  flow  for 
given  depth  ratios  a.  Figure  5  depicts  the  variation  of  neutral  stability  pairs  m.  A  for  three 
different  depth  ratios.  In  each  case  there  is  a  global  maximum  oscillation  amplitude  above 
which  the  flow  is  linearly  stable  for  all  viscosity  ratios.  For  example  when  the  depth  ratio  is 
a  =  1.4  we  see  that  the  flow  is  stable  if  A  >  0.6412  (approximately). 

Finally  we  quantify  the  effect  of  varying  the  frequency  of  the  oscillations,  through  the 
non-dimensional  parameter  fl.  Increasing  the  frequency  reduces  the  effect  of  the  unsteady 
terms  and  the  stability  of  this  fluid  regime  becomes  comparable  with  that  of  steady  plane 
Couette  flow,  (denoted  by  the  broken  line).  These  stability  results  ran  be  understood  by 
consideration  of  the  unperturbed  flow  at  high  oscillation  frequencies.  In  such  a  regime  the 
flow  separates  into  a  Stokes  layer  in  the  vicinity  of  the  oscillating  wall,  and  away  from 
this  layer  the  flow  is  steady  and  corresponds  to  two  phase  (louette  flow  due  to  a  boundary 
which  moves  with  constant  velocity.  As  long  as  the  Stokes  layer  is  thin  compared  to  the 
thickness  of  the  upper  layer,  therefore,  the  interfacial  mode  is  expected  to  be  insensitive 
to  the  wall  modulations  as  indicated  by  our  numerical  results.  The  Stokes  layer  thickness, 
in  nondimensional  terms,  is  proportional  to  and  for  the  highest  values  of  i'l 

used  in  Figures  6(a)-(b)  the  ratio  between  the  Stokes  layer  thickness  and  the  distance  of 
the  unperturbed  interface  from  the  wall  is  approximately  0.025  and  0.25  r(’S])eclively  for  a 
unit  Reynolds  number,  confirming  the  arguments  given  above.  On  the  other  hand,  as  O  is 
reduced  the  stal)ilizing  or  destabilizing  tendencies  characterized  hy  A  >  0.  are  <'mphasiz('d. 
as  sliown  by  Figures  6(a)  and  6(b).  It  can  lie  concluded,  therefore,  that  frum  a  practii-al 
point  of  view,  staliilization  due  to  time  harmonic  modulations  is  likely  to  occur  as  long  as 
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the  modulation  frequency  is  not  too  large. 


5  Conclusions 

We  have  considered  the  effect  of  the  inclusion  of  time  harmoni<’  modulation.s  in  the  un|jei- 
turbed  velocity  field  of  two  phase  (’ouette  How  of  different  liquids.  Our  main  aim  is  the 
demonstration  of  the  stabilizing  effect  that  such  modulations  can  have  on  otherwise  interfa- 
cially  unstable  flow's,  a  finding  that  can  have  useful  practical  appfication.s.  To  tijis  eiid.  we 
considered  the  unsteady  partial  differential  stability  system  in  the  limit  of  U>ng  wavelength 
perturbations  which  can  be  solved  exactly  by  use  of  Floquet  and  perturbation  theory  to 
yield  analytical  expressions  for  the  first  non-zero  Floquet  exponent  which  goxerns  stability 
or  instability.  The  main  conclusion  of  this  study  is  that  inclusion  of  modulations  with  an 
amplitude  larger  than  a  certain  parameter-dependent  threshold,  can  completely  stabilize 
flows  which  are  unstable  in  the  absence  of  modulations;  for  instance  flows  with  the  more 
viscous  fluid  occupying  a  thin  layer  and  bounded  by  the  moving  w'all  can  be  completely 
stabilized  for  long  waves,  at  least.  At  the  same  time,  modulations  can  produce  an  adverse 
eflect  on  the  intei facial  mode.  Flow  arrangements  which  are  stable  or  unstable  can  become 
unstable  or  more  unstable  respectively.  Finally,  we  emphasize  that  the  conclusions  of  this 
work  are  valid  for  long  wavelength  perturbations  which  are  useful  in  providing  analytical 
solutions  to  a  problem  that  needs  to  be  addressed  numerically  in  general.  Further  more,  the 
physical  conclusions  reached  for  the  long  wavelength  limit  are  expected  to  be  indicative  of 
the  behavior  of  general  wavelength  time  periodic  perturbations. 


Appendix 

For  the  solution  of  the  O  (1)  system  (13a-b),  we  impose  boundary  and  interface  conditions 
given  by  equations  (14a-d)  and  (15a-d).  These  yield  the  following  twelve  equations,  wliieh 
determine  the  constants  Aiq,  A20, ....  Fio>  f'io- 


A\q  —  niA'io  ~  0. 

374io  (1  —  ~  3m/42o  d"  B\q  —  iuB-iq  =  0, 

^10  (1  ~  ~  ^20  +  B\q  (1  —  aY  ~~  B20  =  0. 

3^10(1  —  (i)^  —  iAio  +  —  Cl)  —  2B20  H - ; — r  =  0, 

Bio  —  E-io  =  0, 

F\q  —  FiQ  =  0, 

F\q  +  Dio  =  0, 

A 

Cio  sinh  (/3a)  +  D\q  cosh  (/3a)  +  £10^  +  Fio  —  y  =  0, 

=  ,/3C',ocosh(/3a)  +  ,^D,osinh(/3a)  + £,o, 

,t3m“2  C20  +  £^20  =  0rn~iAK, 

C'iosinh(/3)  +  Diocosh  (/3)  =  C72osinh  +  Z32ocosh  (/3m“2  j  , 

Ciocosh  (/3)  +  Dtosiiih  {0)  =  tn“2  [C2ocosh  +  £)2osinh  {^/37tr2)j  . 

At  0  (a)  the  solution  of  (I8a-b),  requires  a  further  eight  equations,  to  determine  the  constants 
All,  A21, ...,  D^i,  D21' 


A\\  —  771^21  — 

All  ~  1)^  +  (fl  ~  1)^  +  C\\  (a  —  1)  +  Dll  +  Xn  (®)  — 


dx 


3A11  (a  -  1)^  +  2Bii  (a  —  1)  +  Cii  +  (°)  ~ 

— A21  +  D21  —  Cix  +  D21  +  X21  (0)  = 

5X21 


3A21  —  2D21  C2\  + 


dy 


(0) 


Dll  —  D21  4-  Xn  (^)  —  X21  (1)  ~ 


<^11 


dy 

2D„ 


dy 


m 


dUx 


dy  J 


(*) 


o  n  .  5^X11  5^X21,,, 

2mB2x  +  -^-^(1) 


0 

0, 

0, 

0 

0, 

0, 

0, 


=  0. 


dy^  dy^ 

Where  xn  and  X21  are  the  particular  solutions  defined  by  equations  (19a-b). 
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Fig.  1(b)  Steady  Growth  With  Deeper  Lower  Layer 
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